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Abstract 

The exactly solvable model of a one dimensional isotropic XY spin chain is employed to study the thermody- 
namics of open systems. For this purpose the chain is subdivided into two parts, one part is considered as the 
system while the rest as the environment or bath. The equilibrium properties of the system display several 
anomalous aspects such as negative entropies, negative specific heat, negative susceptibilities in dependence 
of temperature and coupling strength between system and bath. The statistical mechanics of this system is 
studied in terms of a reduced density matrix. At zero temperature and for a certain parameter values we 
observe a change of the ground state, a situation akin to a quantum phase transition. 
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1. Introduction 

Ordinary thermodynamics describing the macro- 
scopic phenomenology of homogeneous equilibrium 
systems leads to the following stability conditions: 

dP/dV < (1) 
C v > (2) 

These conditions ensure that the pressure P of a 
system decreases as its volume V increases and that 
the system warms up when absorbing energy at 
fixed volume (positive specific heat at constant vol- 
ume Cy). One of the major achievements of math- 
ematical physics was to show that ordinary matter, 
composed of a large number of electrons and pro- 
tons, behave according toEqs. (HHD [J. 

However, it has recently turned out that systems 
not satisfying the condition @, actually exist in 
nature, and that a corresponding thermodynamic 
description should yet be possible in certain cases 
[H, Q . The prototypical examples of such unstable 
systems are stars which are known to expand and 
cool down as their energy increases 0j- 

Within the canonical ensemble the specific heat 
is related to the fluctuations of energy via the re- 
lation ksCy = T~ 2 (SE 2 ). From this expression 
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it is evident that specific heat is necessarily pos- 
itive in the canonical ensemble. Negative specific 
heat, however may appear within the microcanon- 
ical ensemble. This can happen due to different 
mechanisms: 

• The ergodic properties of the system may de- 
pend on energy. Hence, at different energies, 
different parts of phase space are accessible. 
This may lead to negative specific heat and 
other thermodynamical anomalies. An exam- 
ple was given by Hertel and Thirring [f|. 

• The system might be far from the thermody- 
namic limit. Systems that are stable in the 
thermodynamic limit, such as Lennard- Jones 
gases, may display negative specific if only 
their size is small enough 0, @] ■ 

• Long ranged forces might prevent the thermo- 
dynamic limit to exist at all. But even when an 
equilibrium state exists such systems remain 
nonextensive and may show negative specific 
heat [1,0. 

In this work we will present yet another mecha- 
nism that leads to negative specific heat and other 

October U, 2010 



thermodynamic anomalies for systems in contact 
with a heat bath. Apart from extensivity and short 
ranged interactions, another assumption is custom- 
arily made in the statistical mechanics of canoni- 
cal systems, namely weak coupling between system 
and its environment. When the coupling to the en- 
vironment is not negligible, violations of condition 
([2]) may appear even if system and environment, as 
a whole, are in a canonical state Hence, yet 

another item can be added to our list of exceptions: 

• Systems that strongly interact with their envi- 
ronment may display thermodynamic anoma- 
lies 

For example, a single free particle, which would 
display a positive specific heat when weakly coupled 
to a bath, may display negative specific heat when 
in strong interaction with the bath [2] . Similar ef- 
fects were observed for a two level system coupled 
to a harmonic oscillator [lOj. 

We assume that the total system S + B, is in 
weak contact with a super-bath that provides the 
temperature concept (T = ks is Boltz- 

mann constant.) This total system is therefore 
described by the canonical statistics e~^ Htot /Z to t, 
where H to t is the Hamiltonian of the total system 
S + B and Z tot = r £re~ l3Htot . Accordingly, its spe- 
cific heat cannot be negative. However, if the cou- 
pling between S and B is non-negligible, the overall 
canonical state does not factorize into the product 
of two canonical states for S and B respectively: 
e-PBft/Ztot ^ e- fiHs e' fiHB /{Z S Z B ). Although 
the system is in perfect thermal equilibrium with 
its environment, it is not in a canonical state and 
therefore negative entropies and violations of the in- 
equalities (|H2j) may occur. These violations though 
are not a sign of any instability because as a part 
of a stable total system, the system of interest itself 
is also stable. 

There are few available exact solutions of open 
systems displaying this kind of anomalies. One 
known example is the damped free particle 0, 0, [Tljj 
and a two-level fluctuator in contact with a single 
oscillator [10]. Numerical investigations of thermo- 
dynamic anomalies in the context of the Casimir 
effect, the multichannel Kondo effect and of meso- 
scopic superconductors containing magnetic impu- 
rities have been reported recently in the literature 
12l Til 14 1 . In this work we study the thermody- 
namics of a subchain of a longer chain of spins in- 
teracting according to the isotropic XY model with 



S, N s 



3,N B 



Figure 1: (Color online) Schematic representation of the 
model studied. A spin chain is immersed in a thermal envi- 
ronment (referred to as the super-bath in the text) at tem- 
perature T = l/(kg/3). The chain is composed of two parts: 
The system of interest S, made of the first left-most N$ spins 
(green squares), and the "bath" B, made of the remaining 
Ng spins (yellow circles) . When the interaction between the 
two subchains is non negligible, anomalies may occur in the 
thermodynamics of the system of interest S. 



HE^. The spin chain is composed of 



free ends 

two parts. One subchain is defining the system S, 
while the rest of the chain comprises the bath, B 
(see Fig. [JJ. Being this an exactly solvable model, 
we are able to analytically find the relevant ther- 
modynamical functions for the subchain. 

Far from being a purely academic problem, the 
study of the equilibrium properties of spin-1/2 
chains has been recently attracting a great deal of 
attention. Spin systems not only are the basis of the 
physics of magnetic materials [17| but they might 
have an enormous impact with regard to the devel- 
opment of quantum technologies 0, US HH, 22 1- 
We believe that the study of the thermodynamics 
of small quantum systems will help to understand 
Pro blems of quantum information and vice versa 

In Sec. [5] we briefly review the generalities of 
the thermodynamics of open systems. The specific 
model studied in this paper is described in Sec. 121 
while its thermodynamics is illustrated in Sec. 21 
Various anomalies ranging from negative entropy 
and negative specific heat, to negative susceptibil- 
ity are observed for system size, Ns, of the order 
of unity. In Sec. [5] we study the reduced density 
matrix of the open system S and show that it de- 
parts from the canonical form in the strong coupling 
regime. The spectrum of the reduced density ma- 
trix is analyzed and it is observed that, regardless 
of the strong coupling, a quantum phase transition 
occurs at zero temperature. Conclusions will be 
drawn in Sec. El 



2. Thermodynamics of open systems 

Consider the following Hamiltonian: 

Htot — Hs + Hb + Hsb (3) 



2 



describing a systems of interest S interacting with 
a second system, the bath B, via the interaction 
energy term H$b- The compound system S + B 
is in weak contact with a large super-bath at tem- 
perature T = l/fcg/3. Hence it is described by the 
canonical density matrix: 



Ptot 



,-PH t „t 



/Zt, 



where 



Tre 



(4) 



(5) 



is the total system partition function, with Tr de- 
noting the trace over the total system Hilbert space. 
The partition function of the open quantum system, 
S, is given by the ratio of the total system partition 
function Z tot and the bare bath partition function 
Z B , i.e. @, i, M, M, Ezl 0, Q: 



Z — Ztot/Zi 



where 



Tr B e 



(6) 



(7) 



with Tr# denoting the trace over the bath Hilbert 
space. According to the rules of statistical mechan- 
ics, the Helmholtz free energy of the open quantum 
system is then: 



F = -p- x \nZ . 



(8) 



From this free energy, all the relevant thermody- 
namic functions of the open quantum system can 
be derived. In particular, the entropy and the spe- 
cific heat are given by S = —dF/dT, and Cy = 
TdS/dT, respectively. Thanks to Eq. © the sys- 
tem's free energies is the difference of the free en- 
ergy of the total system and the bare bath, i.e.: 



F = F t , 



F E 



(9) 



Due to the linearity of the derivative, this implies 
according relations for the entropy and the specific 
heat of the system, which take the form: 



S — Stot 
C = Ctot 



Sb , 
C B 



(10) 
(11) 



While Ctot and Cb are positive numbers, there is 
no reason why their difference should be positive as 
well. Hence, coupling the bath B strongly to the 
system S may result in an overall decrease of spe- 
cific heat. The same holds for the entropy and other 
thermodynamic quantities that linearly depend on 
the free energy. These relations will be exemplified 
with the model described below. 



3. The Model 

We consider the isotropic XY model of N inter- 
acting spins on a one dimensional lattice of equally 
spaced sites with free ends: 

h N T 
^=2E^ + ^EK^+W (12) 

where erj , and cj are the Pauli matrices of the 
j th spin. Dropping the label j these are: 



1 

1 

-i 

1 

1 

-1 



(13) 
(14) 
(15) 



The first term in the Hamiltonian ([12)) accounts for 
the Zeeman interaction between each spin and an 
applied magnetic field h pointing in the z-direction. 
The second term describes the nearest neighbor in- 
teraction with strength J. 

The isotropic XY model (fl"2"|) is an exactly solv- 
able model. By means of the Jordan- Wigner trans- 
formation [lfjj], the Hamiltonian (fl"2"|) is mapped 
onto a system of N free fermionic eigenmodes, of 
energies 



\{ N) =h-2Jcos, 

1 N + 1 



kn 



1...N, (16) 



see |Appcndix A Thereby, the spectrum of the 
Hamiltonian (fTS"]) is expressed in terms of the oc- 
cupation numbers of each fermionic mode Tlj , &-S 



N 

--n!,...n N = y ] A 



Nh 



(17) 



where each rij can only be or 1. Given 
the spectrum, the partition function, Z^ = 
YL „ e -fei,...n K j s calculated as: 



A' 



z 



N 



-PNh/2 



n > 



fc=i 



e-^H. (18) 



According to the formula F = —f3~ 1 \nZ, the free 
energy of an isotropic XY chain of N spins then 
becomes: 



(N) 
k 



(19) 



fc=i 



3 



where 



AN) 

Jk 



-r 1 ln( 



h 
2 



(20) 



denotes the contribution from the k th fermionic 
mode. 

3.1. A part of the chain as an open system 

We now consider the part of the chain consisting 
of the first N$ spins, counted, say, from the left end 
of the chain, as the system of interest, S, and the 
rest of N — Ns spins as the bath, B, see Fig. Q] The 
total Hamiltonian H^, can be recast in the usual 
system-bath form Hn = Hs + Hb + Hsb of Eq. 
©, with 



N s 



N s -1 

E 



• (21) 



H„ = - 



N S 



j=N s +l 



N-l 

E 

j=N S +l 



( „x x 



Hsi 



J 



\ N> 



N s a N s + l 



"Ns a N s - 



l) 



(22) 
(23) 



The regime of strong coupling holds when the sys- 
tem energy is comparable to the interaction energy. 
Weak coupling is achieved in this problem when at 
least one of the following holds Ns ~> 1, h 3> J, 
ksT 3> max(/i, J). 

From the Eq. @, the thermodynamic free energy 
of the open system S is the difference between the 
free energy of the total chain of length N and the 
free energy of the bare bath £>, which is itself a 
XY chain of a certain length Nb- Its free energy 
is then readily obtained by replacing N with Nb in 
Eq. (flT))) . We thus obtain for the system of interest 
S: 

F = F n -F Nb . (24) 
Using Eqs. ([r3l2"Ul) : 



F 



N 

fe=i 



N B 

E 

9=1 



Jq 



(Nb) 



(25) 



4. Thermodynamics of the subchain S 

We discuss now the explicit results for the rele- 
vant thermodynamical functions, such as entropy, 
specific heat, magnetization and susceptibility of 
the open system defined in Eqs. Q and (|2"T1) . See 
also Fig. [TJ 
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Figure 2: (Color online) Entropy per spin for Ng = 2, TVe = 
18 (black solid line), N s = 18, N B = 2 (red solid line), 
N s = 2, N B = (black dashed line), iV s = 18, N B =0 (red 
dashed line). Top panel: h = 3 J/2, bottom panel h = J/2. 
Entropy becomes negative at low temperature for the lower 
magnetic field h = J/2, in the strongly coupled case, i.e. 
N s = 2, Nb = 18, h = J/2 



4-.1. Entropy 

We begin our discussion with studying the en- 
tropy of iS: 



S= - 



dF 
df 



From Eq. (|19|) . follows 



•>jv 



N N 
(N) 



E' 

fc=i 



E* 

9=1 



(JVb) 
9 



(26) 



(27) 



where SV and Sn b denote the total system and the 
bare bath entropies, respectively, and 



S W = ln(l + e^ ) )+_ 



(N) 



+ e 



0K 



(28) 



is the single mode entropy of the total chain of 
length N. Likewise s q Ns ^ is the single mode en- 
tropy of the bath B. In Fig. [2] the temperature 
dependence of the entropy per spin S/N$ of the 



4 



system S, is shown for different system sizes Ns 
of S and fixed total size N. The magnetic field is 
h = 3 J/2 and h = J/2 in the top and bottom panel 
of Fig. [21 respectively. For comparison, the graphs 
showing the temperature dependence of entropy for 
the same system sizes in absence of the bath, are 
shown as well. We shall use the label So for the 
bare system. 

Notably the presence of the bath does not affect 
much the entropy of the system at high tempera- 
tures, whereas its effect is most prominent at low 
temperatures. This high temperature behavior can 
be understood by looking at Eq. (|2"T)) . The entropy 
sjj. of each mode tends to In 2 for large temper- 
atures, thus the entropy of the total system tends 
to A^ln 2. Likewise, the entropy of the bare bath 
tends to Nb In 2 and that of the system S tends to 
(N — TV's) In 2, which is the same as the entropy 
of the bare system Sq at high temperature, i.e., 
7V,s In 2. Thus, for increasing temperature the ef- 
fect of the bath on the entropy of the system S 
becomes less relevant. 

With reference to both panels of Fig. [5J we see 
that at low temperature, the effect of the bath be- 
comes increasingly relevant as the system size de- 
creases. As the system shrinks, the relative effect 
of the interaction with the bath, which is a surface 
effect, becomes more important (strong coupling). 
In particular, we notice the pronounced peak in the 
entropy for the smallest value of Ns (Ns — 2) in 
the case h = 3 J/2, top panel of Fig. [2] A different 
situation occurs in the bottom panel of Fig. [5]where 
the magnetic field is chosen to be h = J/2 and all 
other parameters are kept unchanged as compared 
to Fig. [21 The major difference is the presence of a 
region of negative entropy. This is a feature that is 
often found in open systems, in the regime of strong 
coupling [li, [11 [13 

Fig. [3] shows the entropy of a larger chain (Ns = 
100) coupled to an equally large bath (Nb = 100) as 
a function of coupling strength J/(ksT), as com- 
pared to the entropy of the same chain with no 
bath. For small enough J/(ksT) the two entropies 
are approximately equal, while thir difference is ap- 
parent for large enough coupling. The entropy of 
the uncoupled chain is always positive, whereas the 
entropy of the coupled chain may become negative 
for large enough coupling J/(ksT). 

4-2. Specific heat 

The specific heat at constant magnetic field h 
is obtained from the entropy via the standard for- 
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Figure 3: (Color online) Entropy per spin as a function of 
J/(kgT) for chain of Ns = 100 spins coupled to a bath of 
Ng = 100 spins (black solid line), as compared to an isolated 
chain with Ns = 100 (red dashed line) at fixed h = 5kgT. 
The two entropies coincide only for small enough values of 
J/(kgT). While the entropy of the uncoupled chain is al- 
ways positive, the entropy of the coupled chain may become 
negative for large enough values of J/(kgT). 



mula: 



C ~ 1 dT 



(29) 



Evidently C is the difference between the total sys- 
tem specific heat, CV, and the bare bath specific 
heat, Cn b 



C = C 



N 



c 



N 
k=l 



,(N) 



N B 

q=l 



(30) 



where the single mode specific heat is: 



JN) _ rp ggfc 



mshGSAfV^ 



(31) 



The specific heat of each single mode vanishes at 
high temperature and so do the total chain specific 
heat, the specific heat of the bath, the specific heat 
of the system S, and the specific heat of the bare 
system Sq. 

Fig. @] shows the specific heat for the same pa- 
rameters as in Fig. H that is N = 20, A^ 5 = 2, 18, 
with h = 3 J/2 (top panel) and J/2 (bottom panel). 
Quite interestingly, for the smallest system (A^s = 
2) which is more affected by the bath, the specific 
heat displays a very pronounced positive peak at 
very low temperature followed by a negative peak at 
intermediate temperatures, for the case h — 3J/2. 
This appearance of a negative specific heat at in- 
termediate temperature was observed also for a free 
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Figure 4: (Color online) Specific heat per spin for Ng = 
2, Ng = 18 (black solid line), N s = 18, Ng = 2 (red solid 
line), N s = 2,N B = (black dashed line), N s = 18, N B = 
(red dashed line). Top panel: h = 3 J/2, bottom panel 
h = J/2. Negative specific heat regions appear at low tem- 
perature in the case h = J/2 (bottom) and at intermediate 
temperatures in the case h = 3 J/2. 



particle coupled to a minimal bath composed of a 
single oscillator [2j. 

The behavior changes if the magnetic field is de- 
creased to a value smaller than J, e.g. h = J/2, 
see bottom panel in Fig. @] In case of very small 
Ns, i.e. Ns = 2, there is a negative peak at very 
low temperature and a positive one at intermediate 
temperature. A similar situation was also observed 
for the case of a two-level system coupled to a min- 
imal bath composed of a single oscillator [Io| . 

The sign of the specific heat is given by the sign 
of the derivative of S with respect to temperature 
T . Thus, whenever the entropy displays a region 
where it decreases for increasing temperature, cor- 
respondingly the specific heat displays a region of 
negative values. Therefore, the presence of a posi- 
tive peak in the entropy at low temperature leads 
to the appearance of negative and positive peaks in 
the specific heat at low and intermediate T's, re- 
spectively (compare the top panels in Figs. [5J UJ. 
Likewise the presence of a negative peak in the en- 



tropy at low T leads to a negative specific heat re- 
gion at low temperature whose extension is smaller 
than that of negative entropy (see also [l0|). 



4-3. Magnetization 

Of further interest is the magnetization of the 
system S as a function of temperature T and mag- 
netic field h: 

dF 
~dh 



M = — 



(32) 



It is given by the difference of the magnetizations 
of the total chain and the part of the chain repre- 
senting the bath, reading: 



M = M N - M Nb = 



N N 

W _\ 

k Z^""9 

fe=l 5=1 



77?., 



??? 



(Nb) 



, (33) 



where the single mode magnetization, rojj, i s: 

mf^tanh^f^) . (34 ) 

At high temperature the single mode magnetiza- 
tion tends to zero as the thermal agitation wins over 
magnetic ordering. Fig. [5] shows the magnetization 
per spin for the same parameters reported in Figs. 
II H i.e., N = 20, N s = 2, 18 and h = 3J/2 (top 
panel), h — J/2 (bottom panel). 

With reference to the top panel, we see that in 
absence of a bath the magnetization per spin of 
the shortest chain, Ns = 2, is much larger than 
the magnetization of the longer chain Ns = 18, for 
thermal energies (fcsT) up to the order of J. When 
the system is put in contact with the bath, no rele- 
vant change in the magnetization is observed for the 
longer chain Ns = 18. For the smallest chain how- 
ever, a quick drop of the magnetization occurs with 
increasing temperature, and values close to those 
pertaining to larger chains are reached already at 
temperatures of the order J/(10fcs)- 

The bottom panel of Fig. [S] shows the magneti- 
zation per spin for the same parameters as of the 
top panel but for h = J/2. In contrast to the case 
h = 3J/2, and with reference to the smallest chain, 
Ns = 2, we observe a dramatic enhancement of sev- 
eral orders of magnitude of magnetization at values 
of ksT/J < 0.1, due to the coupling to the bath. 
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Figure 5: (Color online) Magnetization per spin for JVg = 
2,N B = 18 (black solid line), N s = 18, N B = 2 (red solid 
line), N s = 2,N B = (black dashed line), N s = 18, N B = 
(red dashed line). Top panel: h = 3J/2, bottom panel h = 
J/2. 
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Figure 6: (Color online) Susceptibility per spin for Ng = 
2,N B = 18 (black solid line), N s = 18, N B = 2 (red solid 
line), N s = 2, N B = (black dashed line), N s = 18, N B = 
(red dashed line). Top panel: h = 3J/2, bottom panel h = 
J/2. 



4-4- Susceptibility 

By taking the partial derivative of the magnetiza- 
tion M with respect to magnetic field, one obtains 
the magnetic susceptibility: 



X = 



dM 
~dh 



One finds: 



JV 



Ne 



x = x N -x NB =Y J xi N) -Y,x\ NE 



k=l 



q=l 



where the single mode susceptibility reads: 

P 



(JV) 

Xk 



4cosh 2 (/3A^V2) 



(35) 



(36) 



(37) 



Fig. |5] displays the susceptibility for the same set 
of parameter values used in Figs. 0] [SJ i.e., N = 
20, N s = 2, 1 8 and h = 3 J/2, h = J/2. Note the 
enhancement of susceptibility at low temperature 
in both cases due to the presence of the bath. As 



expected, this is more pronounced for the shortest 
chain. 

Of particular interest is the behavior of the sus- 
ceptibility for values of h < J and the smallest chain 
(bottom panel of Fig. E]). Note the large negative 
enhancement at low temperature for Ns — 2. Not 
only the entropy and the specific heat may display 
anomalous behavior in open systems and take on 
negative values, but also the susceptibility may do 
so. 

4-5. Internal energy 

The system internal energy is obtained from the 
thermodynamic relation: 



U = F + TS 



(38) 



Evidently, also the internal energy of the system is 
given by the difference of the respective quantities 
for the total system and the bare bath, resulting to: 



u = u N - u Nb = y u 



N N B 
(JV) _ \ - 
fc 



(39) 



k=l 
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where the single mode energy is given by: 



,(N) 



-fl N) +Ts^'. (40) 
Using Eqs. (|20I28D for fk and Sk we obtain 
u^=Mi N) \{ N) +h/2 



(AT) 



(41) 



where 



Mr - 



i 



i 



(42) 



is the Fermi distribution giving the average occu- 
pation number of the fc'th fcrmionic eigenmodc. 

In Fig. [7] we show the curve (U(T, h), S(T, h)), 
parameterized by the temperature T, for the fixed 
parameters N = 20, N s = 2 and h = 3 J/2 (top), 
h = J/2 (bottom). Note that the relation of en- 
tropy and internal energy is not one-to one every- 
where, but there is a region of internal energies be- 
longing to different entropies and vice- versa. This is 
a specific feature of strongly coupled open systems, 
which may not appear in the weak coupling regime. 
The graph (U(T,h), S(T,h)) contains many infor- 
mations. Each point on the plotted curves corre- 
sponds to a given temperature T, which is identical 
to the slope of the graph at the very same point. 



dU dU/dT d(F + TS)/dT 



OS dS/dT 



dS/dT 



T 



(43) 



Cusps appear in the curve (U(T,h), S(T,h)) for 
those values of T for which the entropy S(T, h) has 
a local extremum. By comparison with Fig. [2l we 
see that for Ns = 2, h = 3J/2, first the entropy 
goes through a maximum and subsequently a min- 
imum, as the temperature is increased. These two 
extrema correspond to the two cusps one encoun- 
ters in the graph (U(T, h), S(T, h)) as one follows 
the curve in the U-S plane starting at S = 0. Note 
that the curve continuously changes its slope when 
passing through the cusps. 

The bottom panel of Fig. [7] shows the curve 
(U(T,h),S(T,h)) N = 20, N s = 2 and h = J/2. 
Only one cusp appears here in correspondence to 
the single extremum (a minimum) of the graph 
S(T, h), compare to Fig. H 

In presence of vanishingly weak coupling, the spe- 
cific heat is positive, meaning that the entropy is a 
strictly increasing function of T. This precludes the 
possibility of having local extrema in the entropy, 
which in turn excludes the appearance of cusps in 
the S, U plot. 
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Figure 7: Parametric plot of entropy versus internal energy. 
Cusps appear both at h = J/2 (bottom), and h = 3J/2 
(top). For both values of h = J/2 (bottom) and h = 3J/2 
(top) cusps appear where the entropy has a maximum as a 
function of temperature. 



5. The system reduced density matrix 

In this section we study the density matrix of the 
system of interest S. We recall that a super-bath 
thermalizes the total system S + B. This means 
that the total system is in a thermal Gibbs state: 



Pn 



N 



(44) 



The density matrix p of S is obtained by tracing 
out the bath degrees of freedom from p^: 



p = TtbPn ■ 



(45) 



For the sake of simplicity, we limit our discussion to 
the case of a small system with Ns — 2, where the 
effects of the coupling to the bath B are maximal. 
In this case of a system composed of two spins the 
reduced density matrix can be calculated by means 
of two-point correlators, according to the general 
formula 



P = 



(°1 °2 / a l a 2 



(46) 



a, 7— 0,x,y,z 



s 



where er" = 1$, denotes the identity operator of 
the Hilbert space of the i th spin (i — 1,2), and {•) 
denotes quantum expectation values with respect 
to pN given in Eq. (|44[) . In the present case it is 



(47) 



with Tr being the trace over the total system's 
Hilbert space. Using Eq. (|47|) . we find the reduced 
density matrix (in the basis {| + +), | + — ), | — +), | — 
— )} of the common eigenvectors of a\ and erf), as: 



P= 7 



where 
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ai2 
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a 33 ) 



(48) 
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Figure 8: (Color online) The low temperature spectrum of 
the system reduced density matrix consists of eigenvalues 
that change step-like as a function of applied field h/J. 
Parameters are: J/3 = 2000, = 2, N B = 20. Beyond 
\h/J\ = 2 only a single nonvanishing eigenvalue exists. At 
higher temperatures the steps are washed out. 
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a-33 


= 1 


- K) - 


(4) + 
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ai2 


= 2 








(53) 



Using Eq. (|47l) one obtains, after standard but te- 
dious algebra: 



/ „x _x\ 
(°1 °2 ) 



AT- 



-I . / kir \ . ( 2kn 

sin — — - sin 



-iE sir 

k=l 



N+l 



N + lJ k 
(54) 



k—1 



jkir 



(55) 
(56) 



Diagonalization of p is straightforward. We obtain 
the following eigenvectors 



|l)=l++> 

1 2) =sin(9_| 

1 3) = sin 6+ 1 
|4)=|--> 



-) + cos( 
+) + cos( 



where the phases 9± become 

(erf) -((Tf) ±6 



9± = arctan 



2(afaf; 



(57) 
(58) 
(59) 
(60) 



(61) 



The corresponding eigenvalues are: 



Pi 

l>2 
P3 
Pi 



= (1- 

= (H 

=«33 



K^))/4 
K4))/4 



where we have introduced the abbreviation: 



(62) 
(63) 
(64) 
(65) 
(66) 



(67) 
(68) 



5.1. The zero temperature limit 

It is interesting to study the spectrum of the re- 
duced density matrix at zero temperature, i.e., for 
(3 — > oo. Fig. [5] shows the spectrum as a function of 
the applied magnetic field for (3 = 1000, J = 1, and 
N = 20. For h>2J the only nonzero eigenvalue is 
j>4, meaning that the open system S is in the pure 

state 4) = | ). The same holds for h < —2 J too, 

in that case the ground state is |1) = | + +)• This 
is an interesting results: even though the system is 
(strongly) coupled to the bath £>, its ground state 

is the same pure state | } (or | + +)) that the 

system would have in absence of the bath, if the 
magnetic field is strong enough. The fact that the 
bath does not alter the ground state for \h\ > 2 J 
is another interesting aspect of a system strongly 
coupled to its environment. This happens because 
for \h\ > 2 J the total system ground state is a pure 



factorized state given by the product of single spin 
states all pointing in the same direction (parallel to 
h for h < — 2J, antiparallel to h for h > 2J). By 
tracing out the bath spins, a pure factorized state 
given by aligned spins remains for the subsystem 
S. On the other hand, when \h\ < 2J, the spins in 
the ground state of the total system are entangled. 
Thus tracing out the bath degrees of freedom leads 
now to a mixed state for the subsystem S. Note 
that the transition of the system density matrix 
takes place at the same parameter value h/(2J) = 1 
at which the infinite chain (N = oo) undergoes a 
quantum phase transition [30] . 

The abrupt steps appearing in the spectrum, as 
displayed in Fig. [51 correspond to the values of 
h for which the energy of a fermionic eigenmode 
vanishes, i.e, \ { k N) = h - 2Jcos(fc7r/(A + 1)) = 0, 
see Eq. (TIT)]) . These steps stem from the terms 
Afjf^ in Eqs. (|54l55p . which at zero temperature 
tend to unit steps: 

lim Afj; N) = lim 1/(1 + e<°) = 6{-\[ N) ) 

j3—>oo (3— too 

(69) 

where 9 is the Heaviside function. The summation 
over all these steps, then originates the staircase 
structure of the spectrum. In the thermodynamic 
limit Nb — > oo, both width and height of the steps 
shrink and a continuous curve results. An analyti- 
cal expression for the PiS can then easily be found 
by replacing sums with integrals (not reported). 

5.2. The high temperature limit 

It is apparent that, in general, the re- 
duced density matrix is not of the form po = 
e~@ Hs /(Trge~P Hs ), corresponding to the uncou- 
pled case. In the same basis {| + +), | H ), | — 

+), | )} of Eq. (|35|). this canonical p reads: 



Pa = 



1 



( 



-/3h 







V 



with Qq 





cosh(/3 J) - sinh(^J) 
- smh(/3J) cosh(^J) 


(70) 

Tr s e~ 13113 = 2[cosh(/3J) + cosh(/3/i)]. 
The analysis carried out in the previous section 
shows that, at high temperature, the thermody- 
namic behavior of the system is not affected by the 
coupling to the bath. This suggests that at high 
temperature the reduced density matrix p should 
tend to the uncoupled system density matrix pq. 
This is indeed the case. Taylor expansion around 



/3 = reveals that p and po coincide up to second 
order in /3. 

6. Conclusions 

The (possible) appearance of thermodynamic 
anomalies due to a violation of the the usual weak 
coupling assumption have been highlighted. Nega- 
tive specific heats and entropies were reported al- 
ready in the literature 0, 0, [l(j , here, for the first 
time, we described anomalies of the susceptibility. 
Apart from the fact that the susceptibility may be- 
come negative, its value can also be enhanced by 
several orders of magnitude, in the strong coupling 
regime. 

The coupling strength of the XY model is de- 
termined by the exchange energy J. This energy 
has to be compared with thermal energy and the 
bulk energy of the system, which grows with the 
system size Ng. Therefore, the anomalies may 
only emerge for relatively small systems at suffi- 
ciently low tempratures. For large systems and 
at high temperatures, the canonical density matrix 
exp[— Hg/ksT]/Zg describes the equilibrium prop- 
erties of the system excluding the presence of any 
anomaly. This is a well known generic feature of 
any spatially extended system interacting with its 
environment by short range interactions 3l|. 

The analysis of the reduced density matrix shows 
that the thermal equilibrium properties of open 
quantum systems may grossly differ from those re- 
sulting from a canonical Gibbs state. This is true 
even in the thermodynamic limit of the bath Nb — > 
oo, as long as the coupling remains strong, (J /ksT 
is sufficiently large and Ng sufficiently small). In re- 
gard to the case Ng = 2, we found that the ground 
state changes from a pure to a mixed state, a situ- 
ation akin to a quantum phase transition [301 ]. This 
transition takes place at the same value h = 2 J of 
the magnetic field at which the usual phase transi- 
tion in the closed isotropic XY model, is observed. 

Finally we note that for the even smaller sys- 
tem composed of a single spin the above described 
anomalies continue to exist. We have restricted our- 
selves to the properties of spin chains and therefore 
disregarded the case of an open single spin system. 
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Appendix A. Solution of the isotropic XY 
model 

In this appendix we briefly review the solution of 
the isotropic XY model, Hamiltonian (fT2)) [TBI . 16 1. 
New operators are defined by means of the Jordan- 
Wigner transformation (l5j : 



3-1 

n 

fc=i 



(A.l) 



where er. is the j spin lowering operator 



1 



(A.2) 



The operators Oj are fermionic operators satisfying 
the canonical anti-commutation rules: 



{aLak} = 6ik {a,j,a k } = 0. 



(A.3) 



In terms of these operators the Hamiltonian (TT21 is 
expressed as: 



N 



N-l 



Nh 



H N = -h^ a]aj - J Y1 ( a ] a 3+i + a ]+i a 3-) + -7T 



3=1 



3=1 



(A.4) 

Next, following known procedures, yet new opera- 
tors are defined, via the discrete sine Fourier trans- 
form: 



JV 



N ■ 



kin 
N + l 



(A.5) 



The new operators also obey canonical fermionic 
anti-commutation rules: 



{b],h} = kk {bj,b k } = 0> 



(A.6) 



Expressing the Hamiltonian in terms of the fefc's, 
one gets: 



JV 



3=1 



(A,, 



i. e. a free fermionic Hamiltonian, with single mode 
energies , Eq. (fT7i| . Denoting the Fock state as- 
sociated to the operators {b k }k=i..N as \n\, . . . n/v), 
that is: 

b\b k \ni, ...n k . ..n N ) = n k \ni, ...n k . ..n N ) 

(A-8) 

with n k — 0,1, the eigenvalues of the Hamiltonian 
read: 

e ni ,...n N =Y,X k N) n J + ^. (A.9) 

3=1 
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